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FOREWORD 


This  document  has  been  prepared  to  provide  a  set  of  graphical  tools  that  may  aid 
the  orbit  analyst  in  the  selection  of  thrust  parameters  to  be  used  during  the  planning  phase 
of  an  orbit  adjust  design  for  satellites  operating  in  small  eccentricity  orbits.  It  is  also  intended 
to  serve  as  the  formulation  document  for  the  NSW C  orbit  adjust  design  software  program 
ORBADJ.  This  is  a  general  software  program  that  may  be  used  to  generate  thrust  parameters 
for  the  initial  planning  phase  of  an  orbit  adjust  design.  This  report  has  been  reviewed  and 
approved  by  Dr.  R.  J.  Anderle. 
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INTRODUCTION 


The  existence  of  mission-related  constraints  and  flight  profiles  requires  that  periodic 
orbital  maintenance  thrusting  be  performed  on  most  artificial  earth  satellites.  This  is  especially 
true  for  low-altitude  satellites,  since  they  are  continuously  perturbed  by  their  interaction  with 
the  atmosphere.  Operational  orbit  adjust  designs  for  earth  satellites  are  normally  performed  in 
an  iterative  fashion  using  coarse  satellite  ephemerides  for  early  planning  and  updated  precise 

ephemerides  for  the  final  design.  Although  the  desired  orbital  modifications  are  known  from  a 

comparison  of  mission  requirements  with  current  mission  performance,  the  analyst  must  still 
select  the  most  expedient  thrust  parameters  for  the  design.  Initial  estimates  f~>r  these  pa¬ 
rameters  are  often  made  by  using  the  results  obtained  from  orbit  perturbation  equations  and 
are  improved  during  the  iterative  design  procedure  mentioned  above. 

The  following  sections  of  this  report  are  concerned  with  the  development  of  formulations 
and  data  that  can  be  used  to  provide  thrust  parameter  estimates.  Such  information  can  be 
used  as  tools  for  the  applied  theorist  as  well  as  the  operational  analyst.  Considered  first  are 
derivations  of  analytic  expressions  that  relate  changes  in  certain  orbital  parameters  produced 
by  an  applied  impulsive  thrust.  These  relationships  are  then  used  to  provide  some  approximate 
results  that  are  useful  for  satellites  in  small  eccentricity  orbits.  A  general  discussion  of  the 

NSWC  orbit  adjust  design  computer  program  ORBADJ  is  provided  in  the  final  section. 


FORMULATIONS 


The  changes  in  orbital  parameters  caused  by  instantaneous  velocity  impulses  can  be 
determined  from  results  obtained  from  the  Lagrange  planetary  equations.1  These  equations, 
expressed  in  their  Gaussian  component  form,  are  given  by 


a  =  2a2 


(1) 


*P.  M.  Fitzpatrick,  Principles  of  Celestial  Mechanics,  Academic  Press,  Inc.,  New  York,  New  York,  1970. 
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v  here  a,  e,  i,  cj,  and  J2  are  the  usual  Keplerian  elements;  n  is  the  earth’s  gravitational  con¬ 
stant;  0  is  the  true  anomaly;  and  u  is  the  true  argument  of  latitude  defined  as 


u  =  6  +  co  (7) 

The  quantities  r,  v,  and  h  are  the  satellite  radial  distance,  speed,  and  angular  momentum, 
respectively.  These  are  defined  as  follows: 


a  (1  -  e2) 
1  +  e  cos  9 


v  = 


(H) 


h  =  \/#ia  (1  -  e2) 


The  acceleration  components  F(  (i  =  R,  I,  C)  appearing  in  the  planetary  equations  are 
defined  by  the  vector  equation 
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F  =  FR  R  +  Fjl  +  FcC  (11) 

where 

*  r  x  r 

c  =  — r 

|r  x  r| 

4 

A  r 
i  = 

|r| 

and 

AAA 

R  =  I  x  C 

Here  r  and  r  are  the  satellite  position  and  velocity  vectors. 

Verification  that  the  thrust  parameters  associated  with  an  orbit  adjust  design  produce 
the  desired  orbital  parameter  changes  requires  that  the  thrust  profile  be  integrated  over  the 
finite  thrust  interval.  A  very  analytically  tractable  approach,  which  yields  quite  useful 
estimates,  uses  an  impulsive  thrust  model  in  which  the  Dirac  delta  function  describes  the 
thrust  profile.  Specifically,  if  the  thrust  components  defined  by 

Fj  =  AV;  5(t  -  tOA),  (i  =  R,  I,  C)  (15) 

where  5  is  the  Dirac  delta  function  and  tOA  is  time  at  which  the  thrust  impulse  is  applied, 
are  substituted  into  Equations  1  through  6  and  the  equations  are  integrated  over  a  small 
time  interval  centered  about  tOA,  the  following  relations  are  obtained: 


(12) 

(13) 

(14) 
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ere  all  quantities  are  preadjust  values  at  tD  A . 


(20) 


(21) 


These  results  can  be  applied  to  those  obtained  from  two-body  analytics  to  obtain  results 
other  parameters  of  interest.  The  changes  in  apogee  distance,  perigee  distance,  orbital 
iod,  and  true  anomaly  are  given  by 


1 


Arp  = 


a  /  1  +  e\ 

AV*  +  2  -  t  —  J  a  +  COS  0)  iV, 


AV'  +  2  y  (fH)  °  '  °0S  6)  AV‘ 


AP  =  3Pa^  ^  AVj 


(22) 


(23) 


(24) 


(1  +  e2)  cos  d  +  2e 
e(l  +  e  cos  0) 


avr 


sin  6 
ev 


AVI 


(25) 


pectively,  where  P  is  the  Keplerian  period 
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A  SUMMARY  OF  APPROXIMATIONS  VALID  FOR 
SATELLITES  IN  SMALL  ECCENTRICITY  ORBITS 


For  sufficiently  small  eccentricities,  the  following  approximations  may  be  made: 


h  «  V/aa 


By  applying  these  approximations  to  Equations  16  through  24,  one  can  form  the  following 
approximate  ratios  between  thrust  components  and  induced  orbital  parameter  changes: 

Aa  P 

—  *  -  (30) 

AVj  ir 


Ae  /  a  \  'h 

\7] 


(-;) 


2  (  -  \  cos  6 


Ai  /a 


ay 


■ft) 


sin  i  ~  /  —  sin  u 
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(35) 

(36) 


AfA  P 

^  -  (1  +  cos  9)  (37) 

AV,  7T 


ArP 

AVr 


P 

—  sin  6 
2  7T 


(38) 


Arp 

AV, 


(1  -  cos  6) 


(39) 


(40) 


ere  P  is  the  period  given  by  Equation  26.  Ratios  for  Aco  and  A 6  are  not  given  because 
the  small  divisor  problem  associated  with  the  small  eccentricity  assumption.  These  ratios 
plotted  for  the  reader’s  convenience  in  Figures  1  through  9  for  several  orbital  periods 
ween  5200  and  7000  sec. 

Based  on  the  results  given  above,  the  following  generalizations  may  be  made: 

•  Only  in-track  thrusting  can  produce  changes  in  a  and  P. 

•  Only  cross-track  thrusting  can  produce  changes  in  i,  ft,  and  u. 

•  Radial  thrusting  produces  equal  and  opposite  changes  in  rA  and  rp. 

•  In-track  thrusting  at  apogee  (perigee)  produces  no  change  in  apogee  (perigee)  distance, 
but  produces  maximum  change  in  perigee  (apogee)  distance. 
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IGNITION  TRUE  ANOMALY  (deg) 

FIGURE  1.  CHANGE  IN  SEMIMAJOR  AXIS  PER  UNIT  IN-TRACK 
VELOCITY  IMPULSE 
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IGNITION  TRUE  ANOMALY  (deg) 


FIGURE  2.  CHANGE  IN  ECCENTRICITY  PER  UNIT  RADIAL 
VELOCITY  IMPULSE 
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C8. 29.13  03/11/63 

PRINT  21#DPCF#CAPC» DVIK#BA 

FORMAT ( /  #  5X#*TG  PRODUCE  A  *#F6.2#*  SEC  PERICC  CHANGE *# /5X  # *AND  A 

1  *#F7.2»*  DEG  PERIGEE  ROTATION  REOUIRES  AN  INTRACK  VELOCITY  INPULS 
2E  OF  *#F7.2#*  METERS  PER  $EC*#/5X#*  AT  A  TRUE  ANOMALY  OF  *#F7.2# 

3*  OEG* ) 

PRINT  22#  J 

FORMAT ( /#5X#*CCNVFRGENCE  OCCURRED  IN  *#I5,*  ITERATIONS*) 

RETURN 

PRINT  31#  J 

FORMAT (///#5X#*C0NVFRGENCE  DIO  NOT  OCCUR  IN  *#I5#*  ITERATIONS* ) 

STOP 

END 

SUPROUTINt  POSVEL  ( A, C E# C I#C L# G> H# 80# PI# R2» P3# VI* V2# V3 #R# VE L  )  0010 

CALL  NWTRFH  ( CL# CE# E# S INE P#C CSE P )  0020 

HSIN-SIN(H)  0030 

HCOS-COS(H)  0040 

GSIN-SIN(G)  0050 

GCCS-COS(G)  0060 

C IS IN*$ IN ( C I )  0070 

CICUS-COS(CI )  0000 

A11-HCOS*GCCS-HSIN*CICOS*GSIN  0090 

A 12— HCOS*GSIN-HSIN*CICOS*GCOS  0100 

A21-HSIN*GCCS+HC0S*CIC0S*GSIN  0110 

A22»HCOS*CICCS*GCOS-HSIN*GSIN  01 ?C 

A31-CISIN*GSIN  0130 

A3  2«C IS IN*GC  CS  0140 

FUN-SQRT(1.-CE*CE)  015C 

Rl-A* (A11*(CCSEP-CE)+A12*<FUN*SINEP) )  0160 

R2«A*(A21*(CCSEP-CE)+A22*(FUN*SINEP))  0170 

R3"A*(A31*(CCSEP-CE)+A32*(FUN*SINEP) )  0180 

R«A*( l.-CE*CCSEP )  0190 

FUNl«(SOPT(eC*A) )/R  0200 

V1»FUN1*< All* (-SINE P)+A12*FUN*C0SEP)  0210 

V2-FUN1*(A21*(-SINEP )+A22*FUN*C0SEP>  0220 

V3-FUN1*(A31*(-SINEP)*A32*FUN*C0SEP)  0230 

VEL»FUN1*S0RT(1.-(CE*CCSEP)**2)  C24C 

RETURN  0250 

END  0260 

SUBROUTINE  NWTRPH  ( 01# B# E 2#S E#C E ) 

NEWTON  RAPHSCN  SOLUTION  TO  KEPLERS  EOUATION  S LBROUT INE 
E1*QLMB*SIMQL) )  /( 1  .-B*COS  (QL  )  ) 

SE-SIN(El) 

C  E «CQS ( El ) 

E2"El*(0L+B*Sfc~El)/(l.”B*CE) 

IF (ABS (E2-E1 )-  l.E-6  )4#4#3 
E  1  -E  2 
GO  TO  2 
S E  "SIM  E 2  ) 

C*«C0S<E2) 

RETURN 

ENO 

/9 
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PRINT  10 

10  FORMAT! ///#1X#*CVI(M/S)*>1 X#+D VC (M/$)*»1X#*DVR(M/S)*#1X#*DA(KH)*#3 

1X»*D£  *  ,3X#*CI(DEG)*,1X>*0W<DEG)*#1X,*D0M(DE6)*#1X»*DU!0EG)*,1X, 

2*0TA(DEG)*»1X#*CRA(KM)*#1X#*DRP(KM  )*#1X»*0P( SEC  )♦» 1X#*TAI< DEG)*) 

PRINT  11>DDVI>DCVC>DDVR.DA»DE>DDI»DDW»DDC#DDU>DDTA»DRA#DRP#DP.TTA 

11  F0RMAT(/#3(lX#F8.3)#lX#F6.2#lX>F7.AilX»F?.4#lX»F7.2>lX»F8.3»lX> 

1  F7.2#1X/F8.3#1X,F6.2#1X,F8.2#1X,F7.2»1X,F8.2> 

RETURN 


20  PRINT  21»DA#DE#DDI#DDW#D00#DDU#DDTA»DRA#DRP»DP#TTA 

21  F0RMAT(/»28X»F6.2»1X»F7.4#1X#F7.4#1X»F7.2*1X»F8.3#1X»F7.2»1X#F8.3> 

1  lX,Ffc.2#lX#F8.2#lX#F7.2»lX,F8.2> 


RETURN 

END 

SUBROUTINE  R  EQC  MDP  CH,  C  AP  C  ) 

•♦♦♦THIS  ROUTINE  COMPUTES  AN  INTRACK  VELOCITY  IMPULSE  REQUIRED  TO 
■♦♦♦PRODUCE  A  DESIRED  PERIOO  CHANGE  AND  PERIGEE  ROTATION. 

C0MM0N/EL£M/A,E*I#W»0 

CONMON/BUPN/TAI#CVI,DVC,DVR 

REAL  I »  NU 

DATA  NU/3986CC.8/#PI/3.1A1 59265359/ 

J-l 

DPC-DAPC*PI/ieO. 

ovc-o.o 

DVR-0.0 
T  A  I  •  0.00 

P  •  2.*PI*( ( (A*A*A)/MU}**0.5) 

10  TA 10  •  TAI 

R-(A*(1.-E*E ) )/ C 1. ♦E*COS  C  TAI >  I 

V2-  HUM <2./R)-(l./A) ) 

V  •  SORT ( V2  ) 

DVI  •  (  OPCF  ♦  ML  )  /  C  3.  ♦  P  ♦  A  ♦  V  ) 

SA  -0 • 5* ( (DPC*E*V)/DVI) 

TAI-ASIN(SA) 

IF (TAI. LT. 0.0)  TAI-TAI42  .*PI 
TST-ABS(TAI-TAIC) 

IF(TST.LT.C.Cl)  GO  TO  20 
J-J41 

IF ( J  *GT .500 )  GO  TO  30 
GO  TO  10 


20  0VIK-1000.*DVI 
8A-TAI*18C./PI 
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V«SQRT (V2 ) 

U-TAI+h 

VP2-(MU/A)*( (l.+E)/(l.-E) ) 

VP»$QRT (VP2 ) 

R  P»A* ( 1  »-E ) 

H* VP*RP 

IF<U.GT.PI2)U-AM0D<U,PI2) 

DRA-( ( K*S INC  TAI )) /V) ♦DVR+2. ♦( A/V )♦(( 1.+E>/<1. -EH *<1. ♦COSCT AI))* 

1  OVI 

DRP»-C ( R* SIN (TAI ) ) / V)*DVR+2.  +  CA /V  I* ( (1 .-E ) / ( 1 ,*E ) I* ( l.-COS ( TAI ) ) * 

1  DVI 

DP*3.*P*A*(V/MU)PDVI 

DA«2.*AtA*(V/PU)*DVI 

DE*((R*SIN(TAI))/(A*V))A0VR*2.*((C0S(TAI)+E)/V)*DVI 
OI-((R*COS(U  )/H)*OVC 

DW  — < ( ( 1 . +E+E ) *COS( TAI )*2.*E )/<£♦< l.+E*COS< TAI »>*V) )*DVR*2.«( SIN (T 
1  AI)/ (E*V>)*DVI-( (R*SIN(U)*COS(I»)/ (H*S IN ( I ) ) ) *DVC 

DON«((R*SIN(U))/(H*SIN(I)))*OVC 

DU  —  (  (R*SIN(II*CQS(I))/(H*SIN(I)))*DVC 

DTA-DU-OW 

RETURN 

END 

SUBROUTINE  CLTD(J) 

•♦♦♦THIS  PROGRAM  OUTPUTS  THE  THRUST  INOUCED  CHANGES  TO  THE  ORBITAL 
‘♦♦♦♦PARAMETERS. 

COMMON/ BURN/TAI,DVI>0 VC, DVR 

COMMON/D£LTA/DRA,DRP,DP,OA,OE,DI,DW,DOM,DU,DTA 
DATA  PI/3.14159265359/ 

DDVI«1C00.*C VI 
DOVC" 1000 •♦DVC 
DO VR*1000 •♦DVR 
DDI  -oi+ieo./Pi 
00VDV*1BC./PI 
DD0«DOM*180./FI 
DDL*Dl'*18o./PI 
DDTA*0TA*18O./PI 
TTA«TAI*180./PI 

IF(J.NE.l)  GO  TO  20 
PRINT  9 

9  FORMAT! ////#  5 X,* THRUST  PARAMETERS  AND  INOUCED  CHANGES  IN  ORBITAL  P 
1 ARAMETERS* ) 


A-8 


-  *  Cm.  fc  k»  r.  >  n,>t.  a  m 
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3  ♦ARGUMENT  OF  PERIGEE  ■♦#G16. 10#*DEGM/5X, 
A  BRIGHT  ASCFNSIQN  OF  ASCENDING  NODE-A# G16. 10#*DEG*/5X» 
5  *7RIE  ANOMALY  •  ♦, G16. 10#*DEG* ) 


RETURN 

END 

SUBROUTINE  MN AN ( TAI# E » AN ) 

♦♦♦♦♦THIS  ROUTINF  COMPUTES  AN  APPROXIMATE  VALUE  FOR  THE  MEAN 
♦♦♦♦♦ANOMALY  ASSOCIATED  WITH  A  GIVEN  TRIE  ANOMALY  TO  THIRD 
♦♦•♦♦3»CER  IN  ECCENTRICITY. 

AN-TAI-2.^EASIN( TAI)*C.75AEPEPSIN<2.PTAI)-0.33333PEPE^EPSIN(3.P 
1  TAD 


RETURN 

ENC 

SUBROUTINE  CUTC<X#Y,Z#XD#Y0#ZD#R*V1 

^♦♦♦♦*THIS  POUT  INF  OUTPUTS  CARTESIAN  POSITION  AND  VELOCITY  COMPONENTS# 
;♦♦*♦♦*$  WELL  AS  THE  ASSOCIATED  RADIUS  AND  SPEED. 

PRINT  10#X#Y#Z*XD,YD#ZD#R#V 
10  FORMATE /» 10X# ♦*  -♦,G16.10»*KMP/10X# 

1  *Y  -+#G16.10,PKMP/10X# 

1  *1  ■♦#  G16.1C»PKM*/10X# 

1  ♦XDUT  •♦#Glfc.lO#PKM  PER  SECV10X, 

1  ♦YCOT  •♦#G16.10#PKM  PER  SEC*/10X# 

1  AZDCT  •♦#Glfc.lO»PKM  PER  SEC*/10X# 

1  ♦RADlUS«*#Glt.l0#*KMP/10X# 

1  ♦SPFFD  ■♦# G16. 10#*KM  PER  SEC*) 

[ 

m 

RETURN 

END 

SUBROUTINE  AC4UST 
C 

I  ♦♦ ***THlS  COiPUTFS  CHANGES  TO  ORBITAL  PARAMETERS  DUE  TO  IMPULSIVE 
[♦♦♦♦♦ VE LOC ITY  CHANGES. 

CGMMON/ELEM/A»£#  1#W#0 
C0NMQN/3l!RN/TAI#CVI#CVC#DVR 

C3MMON/DELTA/DPA#DPP#DP#DA#DE#OI#DW,DCN#DU#DTA 

: 

REAL  I#MU 

C 

DATA  1/3.1  A 159265  359/ 

C 

R*A*(l.-FPE)/(1.4EPC0SITAin 

c 

P«2.*PIA<  t(A#APA)/MU)PP.5) 

PI2*2  »*®I 

V2*Ml*( (2./R  )-( l./A) ) 
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121 


122 


123 

C 

C 


C 

C 

C 

13C 

120 

C 

C 


V-V8R-MJ-1)*VGR 
IF ( IDC .EC .C ) GO  TO  121 
IF( IOC.EQ.l ) GO  TO  122 
DVR-V 
OVI-O.O 
PVC-0.0 
GO  TO  123 
DVI-V 
DVR-0.0 
DVC-0.0 
GO  TO  123 
OVC-V 
DVR-0.0 
DVI-0.0 
CONTINUE 


DO  130  K-ljIJA 

TAI -AB* (K-l) *AG 

CALL  ADJUST 

CALL  OUTO(K) 


CONTINUE 

CONTINUE 


STOP 


C 

C 

C*****ENTER  MODE-2  PROCESSING 

c 

c 

200  PRINT  10 
PRINT  210 

210  FORMATE/  #1X,*M0DE«2  ORBIT  ADJUST  DESIGN*) 
C 
C 

READ  *#DPCH#DAPC 
C 

CALL  REQCH(DPCH,OAPC) 

C 

GO  TO  220 
C 

END 

SUBROUTINE  CUTK ( A,E, INC# PER»RA# ANOM ) 

REAL  INC 


C 

C*****TH IS  ROUTINE  OUTPUTS  KEPLERIAN  ELEMENTS 
C 


PRINT  10#A#E#INC#PER#RA# ANOM 
1C  FORMAT  (/#5X#*SEMIMAJCR  AXIS 

1  *ECCENTR ICITT 

2  INCLINATION 


•*#G16.10#*KM*/5X# 

■*»G16.10/5X> 

■*#G16.10#*DEG*/5X# 
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CALL  ADJUST 
C 

CALL  OUTD(l) 

C 

PRINT  14 

14  FORMAT! // /#5X#*PCST-ADJUST  KEPLERIAN  ELEMENTS  AND  CARTESIAN  VECTOR 
IS*) 

C 

A-A+DA 

E-E+DE 

INC-INC+0I*16O./PI 

I-I+DI 

PER-PER+DW*180./PI 

W-h+DW 

RA-RA*D0M*18C./PI 

0-0*D0M 

ANON- ANOH+DT A* 160. /PI 
TAI-TAI*DTA 
C 

CALL  OUTK (A>  E* INC#  PER.RA. ANOM) 

C 

CALL  MNAN(TA1»E* AH) 

C 

CALL  POSVEL(A*E#I,AM,V,0#MU#X*Y#Z»XD»YD#ZD»R>V) 

C 

CALL  OUTC(X»Y»Z*XD>YD*ZG»R»V) 

C 

STOP 

C 

C 

C*****BEGIN  MODE-1  PROCESSING 
C 

c 

100  READ  *#ANB#ANE«ANG#IDC*VBfVE>VG 
C 

C*****CONVERT  TO  RADIAN  MEASURE  ANO  KM.  PER  SEC. 

C 

AB-ANB*PI/18C. 

AE-ANE*PI/180. 

AG-ANG*PI/18C. 

C 

VBR-VB/1000. 

VER-VE / 1000. 

VGR- VG/10CC . 

C 

C 

PRINT  10 
PRINT  110 

110  FQRNAT(///»1X#* NODE-1  ORBIT  PARAMETER  CHANGE  ENVELOPES*) 

PRINT  10 
C 

IJA-<  (AE-ABJ/AGM2 
I4V-(<VER-VBR)/VGR)*2 
C 

DO  120  4-l»IJV 
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♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦A********************** 


CGMM0N/ELEM/A*E* I*W*0 
COMMON/ BURN/ T A1 *DV I *D VC* DVR 

C  CMMQN/D EL 7A /DRA*DRP* DP* DA*DE*DI*DW*DOM*DU* DTA 
REAL  INC* I*ML 

DATA  P 1/3. 14 1592653 59/* MU/39  8600. 8/ 


READ  ♦  *  MODE 
READ  ♦  *  A* E» INC 
READ  *#  PER* PA 

♦♦♦♦♦CONVERT  DEGREES  TO  RAOlANS 

I-INC*PI/180. 

W-PERAPI/ieC  . 

Q-RA*PI/1P0. 

C 

C*****SELECT  PRCPEP  MODE  LOGIC 
C 

I F ( MODE  .  E 0 • 1 )  GO  TO  IOC 
IE (M0DE.EQ.2  )  GO  TO  200 
C 

(♦♦♦♦♦ENTER  THE  MCDE-0  PROCESSING 
C 

READ  ♦,ANC**CELI*DELC»OELR 
C 

TAI"ANQH* PI/180* 

DVI-DELI/10CC. 

DVC-DELC/100C. 

DVR-DFLR/10CC. 

C 

PRINT  10 

10  FORMAT ( 1H1 ) 

PRINT  11 

11  FORMAT(///*lX**M0DE"0  0R8IT  ADJUST  DESIGN*) 

C 

221  pp ^ 2 

12  F0RMAT(///»5X**PREA0JUST  KEPLERIAN  ELEMENTS  AND  CARTESIAN  VECTOPS* 
1) 

CALL  OUT* ( A#E*INC#PER*RA*ANOM) 

C 

CALL  MNAN(TAI*E*AM) 

C 

CALL  POSVEL(A#E*I*AM»W*C*MU#X*Y»Z*XD*YD*ZD*R*V) 

C 

CALL  OUTC ( X*Y*Z*XD#YD*ZD»R*V) 

C 

C*****COMPUTE  THRUST  INDUCED  CHANGES  TO  ORBITAL  PARAMETERS 
C 
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PRCGRAM  0R8ACJ  ( INPUT#CUTPUT#TAPE6«0UTPUT) 

**•+*«* ************************** *******♦*♦*♦**#**#*#**♦#**#♦***# 
* 

*  THIS  PROGRAM  COMPUTES  THE  CHANGES  PRODUCED  IN  KEPLERIAN 

*  ELEMENTS  AND  OTHER  RELATED  ORBITAL  PARAMETERS  DUE  TO  IMPULSIVE 

*  THRUSTING.  THREE  COMPUTATIONAL  MODES  ARE  AVAILABLE  AND  ARE 

*  DELINEATED  IN  THE  INPUT  CARO  STRUCTURE  BELOH  ■ 


CARD  K  COMPUTATIONAL  MOOE  FLAG  (MODE).  FOR 

MODE-O  ORBITAL  PARAMETER  CHANGES  AND  PRE-  AND  POST¬ 
ORBIT  ADJUST  KEPLERIAN  ELEMENTS  AND  THE 
ASSOCIATED  CARTESIAN  POSITION  AND  VELOCITY 
COMPONENTS  ARE  COMPUTED  BASED  UPON  A  USERS 
DESIGN. 

MODE-1  ORBITAL  PARAMETER  CHANGES  ARE  COMPUTED  FOR  A 
USER  SPECIFIED  DELTA-VELOCITY  RANGE  AND 
GRANULARITY  AND  IGNITION  TRUE  ANOMALY  RANGE 
AND  GRANULARITY. 

MCDE-2  AN  INTRACK  OELTA-VELOCITY  AND  IGNITION  TRUE 
ANOMALY  ARE  COMPUTED  FOR  A  USER  SPECIFIED 
PERIOD  AND  ARGUMENT  OF  PERIGEE  CHANGE.  ALSO 
COMPUTED  ARE  THE  ASSOCIATED  CHANGES  IN  THE 
ORBITAL  PARAMETERS. 

CARD  2t  FIRST  KEPLERIAN  ELEMENT  CARD  CONTAINING  IN  FREE  FORMAT 
THE  SEMI-MAJOR  AXIS  IN  KILOMETERS#  THE  ECCENTRICITY# 

AND  THE  INCLINATION  IN  DEGREES. 

CARD  3t  SECOND  KEPLERIAN  ELEMENT  CARD  CONTAINING  IN  FREE  FORMAT 
THE  ARGUMENT  OF  PERIGEE  AND  RIGHT  ASCENSION  OF  THE 
ASCENDING  NODE  BOTH  EXPRESSED  IN  DEGREES. 

FOR  MOOE-C  ONLY  l 

CARO  91  IGNITION  TRUE  ANOMALY  IN  DEGREES#  IN-TRACK  VELOCITY 
IMPULSE  IN  METERS  PER  SECOND#  CROSS-TRACK  VELOCITY 
IMPULSE  IN  METERS  PER  SECOND#  AND  RADIAL  VELOCITY 
IMPULSE  IN  METERS  PFR  SECOND  • 

FOR  MOOE-1  ONLY  t 

CARD  91  STARTING  IGNITION  TRUE  ANOMALY#  STOPPING  IGNITION  TRUE 
ANOMALY#  IGNITION  TRUE  ANOMALY  GRANULARITY  ALL  IN 
DEGREES)  IMPULSE  DIRECTION  CODE  (-0  FOR  IN-TRACK#>1 
FCR  CR0SS-TRACK#>2  FOR  RADIAL))  ASSOCIATED  STARTING 
IMPULSE#  STOPPING  IMPULSE#  AND  IMPULSE  GRANULARITY 
(ALL  IN  METERS  PER  SECONO). 

FOR  MOOE-2  ONLY  l 

CARD  9t  THE  DESIRED  PERIOD  CHANGE  IN  SECONDS)  THE  DESIRED 
CHANGE  IN  ARGUMENT  OF  PERIGEE  IN  OEGREES. 
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LISTING  OF  PROGRAM  ORBADJ 


impulse  range  and  granularity  and  ignition  true  anomaly  range  and  granularity.  The  last  mode 
computes  the  in-track  velocity  impulse  and  ignition  true  anomaly  required  to  produce  a 
user-supplied  period  change  and  argument  of  perigee  rotation.  Changes  in  the  orbital  parameters 
are  also  computed,  as  well  as  the  associated  preadjust  and  postadjust  Keplerian  elements  and 
Cartesian  position  and  velocity  components.  A  listing  of  this  program  is  provided  in 
Appendix  A.  Samples  of  ORBADJ  output  are  presented  in  Appendix  B. 
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IGNITION  TRUE  ANOMALY  (deg) 

FIGURE  9.  CHANGE  IN  ORBITAL  PERIOD  PER  UNIT  IN  TRACK 
VELOCITY  IMPULSE 

To  illustrate  the  utility  of  Figures  1  through  9,  consider  the  simple  example  where  an 
orbit  adjust  must  be  designed  for  a  low-altitude  satellite  operating  in  a  5200-s  orbit.  It  is 
required  that  this  adjust  increase  the  period  by  40  s  with  no  corresponding  change  in  perigee 
distance.  From  Figure  9  it  is  seen  that  an  in-track  velocity  impulse  of  approximately  20  m/s 
applied  at  any  true  anomaly  will  increase  the  period  by  40  s.  Since  no  change  in  perigee 
distance  is  wanted,  the  thrust  must  be  applied  at  perigee  (i.e.,  a  true  anomaly  of  0°). 
However,  as  can  be  seen  from  Figures  3  and  7,  the  eccentricity  and  apogee  distance  will 
be  increased  by  about  5.1  x  10” 3  and  66  km,  respectively.  Since  cross-track  thrusting  is  not 
applied,  the  inclination,  right  ascension  of  the  ascending  node,  and  true  argument  of  latitude 
will  not  be  changed. 


REMARKS  CONCERNING  PROGRAM  ORBADJ 

The  exact  impulsive  thrust  equations  of  the  Formulations  section  have  been  implemented 
into  an  NSWC  computer  program  called  ORBADJ.  This  program  can  be  used  to  provide 
initial  orbit  adjust  designs  and  can  operate  in  one  of  three  different  user  selectable  modes. 
The  first  mode  computes  orbital  parameter  changes,  pre-  and  post-orbit  adjust  Keplerian  ele¬ 
ments,  and  the  associated  Cartesian  position  and  velocity  components  based  on  a  user-supplied 
design.  The  second  mode  computes  orbital  parameter  changes  based  on  a  user-supplied  velocity 
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—I - 1 - 1 - 

60  120  180 

IGNITION  TRUE  ARGUMENT 

CHANGES  IN  RIGHT  ASCENJ 
AND  TRUE  ARGUMENT  OF 
TRACK  VELOCITY  IMPULSE 


P  =  7000  s 
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